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Abstract 

A new simulation box setup is introduced for the precise description of the wavepacket evolution 
of two electronic systems in intense laser pulses. In this box, the regions of the hydrogen molecule 
H2, and singly and doubly ionized species, and H^^, are well recognized and their time- 
dependent populations are calculated at different laser field intensities. In addition, some new 
regions are introduced and characterized as quasi-double ionization and their time-dependencies 
on the laser field intensity are calculated and analyzed. The adopted simulation box setup is special 
in that it assures proper evaluation of the second ionization. In this study, the dynamics of the 
electrons and nuclei of the hydrogen molecule are separated based on the adiabatic approximation. 
The time-dependent Schrodinger and Newton equations are solved simultaneously for the electrons 
and the nuclei, respectively. Laser pulses of 390 nm wavelength at four different intensities (i.e. 
1 X 10^^, 5 X 10^'^, 1 X 10^^, and 5 x 10^^ W cm~^) are used in these simulations. Details of the central 
II2 region is also presented and discussed. This region is divided into four sub-regions related to 
the ionic state H+H~ and covalent (natural) state HH. The effect of the motion of nuclei on the 
enhanced ionization is discussed. Finally, some different time-dependent properties are calculated 
and their dependencies on the intensity of the laser pulse are studied, and their correlations with 
the populations of different regions are analyzed. 
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I. INTRODUCTION 



Recently, the major focus of the theoretical and experimental research efforts has been put 
on the interaction of two electron systems with ultrashort intense laser pulses pl4T8]. This 
interaction causes many important phenomena such as single and double ionization, charge 
resonance enhanced ionization, dissociative-ionization, and high order harmonic generation 
[in]. Among these phenomena, single ionization has been intensively studied for many years. 
While, double ionization by intense laser fields has continued to draw considerable theoretical 
and experimental attention in the last two decades. The double ionization can occur in 
two ways: i) each electron absorbs energy from the field independently (sequential) or ii) 
one electron absorbs the energy from the field and then shares it with the second electron 
via electron-electron correlation (nonsequential). An important question is about the role 
of electron correlation in the double ionization. Despite its long history, the underlying 
question of the dynamics of electron correlation is still one of the fundamental puzzles in 
quantum physics [H [3, E] . The first step for understanding the behaviour of many-electron 
systems in ultrashort intense laser pulse is to investigate the simplest two-electron systems 
such as the helium atom and the H2 diatomic molecule. One approach to describe this 
interaction is that only one electron is active and responsible for the interaction. This is 
called the single active electron (SAE) approximation. In this approximation, all other 
electrons are assumed to contribute in the dynamics of the active electron through a static 
screening potential. However, for the full description and understanding of the behaviour 
of two-electron systems, it is necessary to consider explicitly both electrons simultaneously. 
This plays a central role in developing our understanding of the interaction of two or many- 
electron atoms and molecules with ultrashort intense laser fields. For this, it is necessary 
to solve the time-dependent Schrodinger equation (TDSE) for the two-electron systems 
that requires a huge amount of computations which is beyond the capabilities of current 
computing facilities. One way to overcome this difficulty is to reduce dimensions of the 
problem, accomponied by the use of soft-core potential [TOl [TTl [231425] . 

In this work, we focus on the electrons dynamics in hydrogen molecule considering the 
nuclei as classical particles. This model reduces the complexity of the problem and helps 
detailing the electrons dynamics without involving the dissociative-ionization process [T3l- 
[T8] . The indistinguishability concept and the symmetry properties of the two electrons will 
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be demonstrated with some details. The paper is organized as follows: In Sec. II, details of 
the simulation box and the numerical solution of the TDSE are described. In Section III, 
results obtained for different intensities of the laser pulses for the introduced simulation box 
are presented and discussed. Finally, the conclusion appears in Sec. IV. We use the atomic 
unit throughout this article unless stated explicitly. 



II. NUMERICAL SOLUTION OF THE TDSE 



Assuming a linearly polarized laser pulse, and considering the fact that dynamics of the 
electrons and the nuclei occur in the laser field direction, we adopt a one-dimensional model 
for both the electrons and nuclei. In what follows, Ri and R2 indicate the nuclei positions 
and Zi and Z2 are the electrons coordinates. Furthermore, M and m indicate the nucleus 
and electron masses, respectively, and e is the electron charge. The temporal evolution of 
electrons of this system is described by the time dependent Schrodinger equation (TDSE), 

i.e. [iniiiiiES] 

.dip{zi,Z2,t;Ri{t),R2{t)) + jD u\ + V i+\ V n^ 

I — =He[Zi, Z2, t; R2{t)}i>{Zi, Z2, t; Rl{t), R2{t)) (1) 

at 

where the electronic Hamiltonian for this system, He(zi,Z2,^; Ri{t) , R2{t)) , is given by 
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^iz,-R^y + aJ y/{z^ - Z2Y + h ^(R,-R,y + c 
+ {zi + Z2)Eof{t) cos(wt), 

(3) 

where Z\ = \ and Z2 = 1 are the charges of nuclei and the screening parameters a, h and c 
are responsible for the softening of the electron-nuclei, electron-electron and nucleus-nucleus 
interactions, respectively. The values of these parameters are set to the same values as used 
by Camiolo et al. [lOj. The laser-molecule interaction energy is formulated in the dipole 
approximation, where Eq is the laser peak amplitude and u = 27rz/ is the angular frequency 
of the laser pulse. The laser pulse envelope, f{t), is set as 

fit) = sin{-7:f (4) 



where ti is the time duration of the field irradiation, set at ri=8 cycles in this work. After 
this time, the simulation continues for 7 more cycles, i.e. with the laser field off, as shown 
in the Fig. [TJ 
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FIG. 1. The laser electric field used in this work has a sine squared envelop function with a duration 
of 8 optical cycles. 

The separated dynamics of the electrons and the nuclei are investigated in quantum and 
classic approaches, respectively. So, we solve the time-dependent Schrodinger and New- 
ton equations simultaneously for the electrons and the nuclei, respectively. This implies 
that adiabatic approximation has been used to separate motions of electrons and nuclei [26] . 
The initial state is a singlet ground electronic state with an equilibrium internuclear distance 
i? = 2.13 at rest. In this singlet state, the fermion electrons adopt an antisymmetric spin con- 
figuration, and thus the electronic spatial part of the wavefunction ip{zi, Z2,t; Ri(t) , R2(t)) 
is symmetric with respect to the permutation of the two electrons. 




III. SIMULATION BOX 

In our previous work, a box is introduced to be used for simulating the behaviour of H2 
in intense laser field [H]. In the present study, we develope a new simulation box shown 
in Fig. [2] Here also, the horizontal and vertical axes of the simulation box are assigned to 
the two zi and Z2 coordinates of the two electrons. Therefore, the ip{zi, Z2, t; Ri(t), R2(t)) is 
symmetric with respect to the principle diagonal (D) of this coordinate system. This sym- 
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FIG. 2. (Color online) The schematic diagram of the simulation box adopted in this study in 
which adjustable and parametrs (with djj+ > ^^2) are introduced. 

metry decreases the processing time and the required amount of memory in the simulations, 
and therefore speeds up the calculations. The simulation is thus carried out only on the 
upper part of the box with respect to this D diagonal. The initial state corresponds to the 
ground state of the hydrogen molecule, i.e. evolution of the system starts from the central 
H2 region. Fig. [2j When an electron of a molecule or ion moves a distance far enough away 
from the nuclei, it becomes ionized. This distance is smaller for H2 (c^/^a) than that for 
{djj+) because the electron of H2 must escape from a core with one positive charge, whereas 
the electron of must escape from a core with two positive charges. To our knowledge, in 
all previous studies, see for example Refs. [IDl EH EH], these two distances {dn-i and c?^+) 
have been assumed to be equal. In this work, we consider them to be different (Fig. |2|. This 
results in the appearance of some new geometries and regions which introduce interesting 
evolution stages for the wavepacket. 

5 



The first ionization of H2 may occur via four pathways, denoted by arrows labled 1-4 
in Fig. [2] As stated in [llj, because of the diagonal symmetry with respect to the D hne, 
pathways 3 and 4 correspond identically to pathways 1 and 2, respectively. In pathway 1, 
electron ei moves away from the nuclei in the negative direction, and in the pathway 2, 
electron 62 gets distance from the nuclei in the positive direction. Therefore, there are two 
distinguishable regions which are labelled by H^(I) and H^(II) (Fig. [2|. 

It is necessary to mention that beyond the H2 region, the behaviour of the ionized electron 
is controlled by the laser electric field. As it is seen in Fig. |2| with introduction of > dnz, 
the regions II^(I) and H^(II) overlap in some spaces, labled as Q(I), Q(II) and Q(III) regions. 
In these overlapping regions, both ei and 62 electrons are ionized. So it seems that these 
regions belong to H^^. However, when ei or 62 belonging to these regions, moves away, it 
enters the corresponding regions. For example, when ei belonging to the upper Q(III) 
region, leaves this region via pathway 5, 62 is captured by the core unexpectedly, and as 
a result, the system enters the upper part of H^(II). This phenomenon can be explained 
as fallows: When the system is in the overlapping Q(I), Q(II) and Q(III) regions, the ei 
and 62 electrons shield each other partially with respect to the charges of nuclei. When ei 
takes distance from the nuclei via pathway 5, its corresponding shielding effect on the 62 
electron is reduced. At this time, the 62 electron which has come now closer to the nuclei, 
introduces more effective shielding on the ei electron so that the ei electron is sought to be 
ionized faster while the 62 electron is pulled closer by the nuclei and thus enters the H^^(II) 
region. Similarly, in the upper Q(III) region, when the ei becomes closer to the nuclei (via 
pathway 7), it introduces a more effective shielding on the charges of the nuclei, and thus, 
the 62 which is now farther from the nuclei, is less affected by the charges of the nuclei. 
Therefore, the 62 will be ionized and thus the ei will be retracted by H^^, and the system 
settles in the upper Hj(I) region. Because of these interesting features of the Q regions, we 
call them as quasi-H^^ regions. The wavefunction in the Q(I) and Q(II) regions are different 
due to the symmetry breaking intriduced by the presence of the linearly polarized laser field. 
Therefore, we have three distinct quasi-H^^ regions. 

In the Q(III) regions, the ei and 62 electrons move away from origin in opposite directions. 
For example, in the upper Q(III) regions, ei is located in the negative part of the zi coor- 
dinate and 62 is located in the positive part of the Z2 coordinate. For the quasi-H^^ species 
in the Q(I) and the Q(II) regions, ei and 62 are located in the same directions with respect 
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to the nuclei, i.e. in the Q(I) and Q(II) regions, ei and 62 are both located respectively in 
the negative and positive parts of the zi and Z2 coordinates. 

Since, all regions of the simulation box are symmetric with respect to the diagonal D, it is 
just necessary to solve the TDSE and calculate the wavefunction for one side of this diagonal; 
the results for the other side can be generated just by switching the ei and 62 coordinates 
in the wavefunction. For example, in the upper H^(I) region, the behaviours of the Ci and 
62 electrons are distinguishable, because ei is ionized and 62 is bound. The same is true for 
lower the H^(I) region, but in which 62 is ionized and ei is bound. However, when both upper 
and lower H^(I) regions are considered together, these distinguishability vanishes. The same 
trends applies to other regions of the simulation box. Therefore, distinguishability of the 
two electrons vanishes if both parts of the simulation box are considered. 

If the intensity of the laser pulse after the first ionization is still strong enough to create 
the second ionization, the species can be ionized and the electrons are introduced to 
the (I-III) regions. In this situation, it is necessary to follow and conserve in the 
simulation box till there is some probability for the second ionization, i.e. the regions 
must be extended so that is not absorbed by the boundary absorption. This condition 
is not met in the setup of the simulation box in Fig. |2] in which the may be absorbed 
by its absorbing boundaries and considered as species, before the laser pulse is turned 
off, and so, some chances of the second ionization is lost. Therefore, the rate of the second 
ionization in the simulation box of Fig. |2] may be underestimated. To cure this problem, 
a new simulation box (Fig. |3| is introduced in which the components is conserved 
until the laser pulse is turned off to assure that the overall second ionization is taken into 
account completely. The regions in this simulation box (Fig. |3| are extended farther 
along both directions of zi and Z2 in comparison with Fig. [2j In this setup of the simulation 
box, when an electron evolves from the regions to the H^^ regions, or it is absorbed 
by the borders labelled 2, 4, 7, and 9, the second ionization takes place. So, the overall 
time- dependent second ionization is equal to the population in the H^^ regions plus the 
population absorbed by the 1, 2, 4-7,9 and 10 absorbing borders in the upper part of the 
box and their corresponding populations in the lower part of the box. Therefore, the overall 
time-dependent population of is equal to the population existing in the regions plus 
the population absorbed by the absorbing borders 3 and 8 and their corresponding borders 
in the lower part of the box. 
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FIG. 3. (Color online) Same as Fig. [2| but the regions are extended to assure that the overall 
second ionization is taken into account. 

We present the results of simulation for two different cases of fixed and freed nuclei in 
the latter of which the dynamics of the nuclei is treated classically. Laser pulses of 390 nm 
wavelength with four different intensities of 1 x 10^^, 5 x 10^^, 1 x 10^^, and 5 x 10^^ W 
cm~^ are used. The shape of the laser pulse is shown in Fig. [ij The time step is set to 6t= 
0.02. The differential operators in Eq. (2) are discretized by the 11-point difference formulas 
which have tenth-order accuracies p9]. To solve the above TDSE numerically, we adopted 
a general nonlinear coordinate transformation for electronic coordinates. For the spatial 
discretization, we have constructed a finite difference scheme with a nonuniform (adaptive) 
grid for the Zi and Z2 electronic coordinates, which are finest near the nuclei and coarsest 
at the border regions of the simulation box. The size of the computational box along both 
Zi and Z2 axes is ±785a.M. and the width of the absorption regions are ±30a.M.. For each 
of Zi and Z2 axes, 4000 points are considered . The interval between these points is 0.2 
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a.u. near the nuclei and becomes ~ 0.75a.u. at the border regions of the simulation box. 
The absorber regions are introduced by using fourth-order optical potentials at the zi and Z2 
boundaries, in order to capture the photoelectrons and prevent the reflection of the outgoing 
wave packets at the borders of the grid. More details of our calculations are described in 
our previous work [30j. Our simulations for different box sizes and grid points show that 
the larger box sizes and more dense grid points does not improve the results. More details 
of our calculations are described in our previous reports [291 EO] • 



SIMULATION RESULTS 



To determine the optimum values of the dnz, d^+, du2^ geometric parameters of the 
simulation box (Fig. |3]), results obtained for simulation boxes with different values of these 
parameters are evaluated and compared. Fig. |4] shows the time-dependent population of the 
H2 region for two different simulation boxes with = 60a. u. once with a box limited (L) 
to the H2 region only with no other regions, and once with the full box (F) as shown in 
Fig. [2j It can be seen from Fig. |4] that variations of the populations for the L and F boxes 
are almost similar. For the comparative study intended here, a fixed size of = 60 for 
the H2 region suffices. However, for precise evaluation of the population transfer between 
regions, this size should be optimized at each intensity. To adjust the and d^2+ values, 
population of the region obtained with different values of these parameters are calculated 
and compared for the cross box (Fig. [3]) and the full box (Fig. |2| and plotted in Fig. |5} This 
figure shows that, at the laser intensity of 5 x 10^^ W cm^^, population of in the cross 
box with d^+ = 500 and dH2_^_ = 580 is converged to that of the full box (with djj+ = 500 
and djj2^ = 1510). The same results are obtained for other laser intensities examined in this 
work. Therefore, for all intensities studied in this work, we set dH2, dfj+ and dH2j^ sizes at 
60, 500 and 580 (i.e. corresponding to the positions of the region boundaries at ±30, ±250 
and ±290), respectively. 



H2 REGIONS 

As shown in Fig. |6| the H2 region can be divided into four sub- regions. The H2(I) and 
(II) regions are related to the ionic state H+H^, while the H2(III) regions are related to 
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FIG. 4. (Color online) Time-dependent population of the H2 region for the limited (L) and full 
(F) simulation boxes obtained at different laser pulse intensities. 
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FIG. 5. (Color online) Time-dependent population of the H2 region for different simulation boxes, 
respectively with djj+ = 200 {d^2^ = 280 (cross) and 1510 (full)), = 300 {dfj2^ = 380 (cross) 
and 1510 (fuh)), d^+ = 400 (d^2_^ = 480 (cross) and 1510 (full)), and = 500 {d^?, = 580 
(cross) and 1510 (full)), at laser intensity of 5 x 10-*^^ W cm~^. 

the covalent (natural) state HH. The time-dependent populations of these different regions 
obtained at different laser pulse intensities for both cases of fixed and freed nuclei are cal- 
culated and plotted in the Fig. [7] . It can be seen from this figure that, at the intensity of 
1 X 10^"^ iycm~^, the ionization is low and the outgoing population from the H2 region is 
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negligible. 




FIG. 6. (Color online) Details of the H2 region of the simulation box introduced in Fig. [Sj 

Under the irradiation of the laser field, the initial symmetric relations between the H2(I) 
and H2(n) regions disappears and, as can be seen in Fig. [7| their electron populations vary 
independently. In the midway of the passage of the exchange of population between the 
H^H+(I) and H+H~(II) ionic regions, populations of the covalent HH(III) regions increase. 
This feature of the wavepacket evolution appears as a small increase in the population of the 
H2(III) region in the time intervals that the populations of the two ionic regions cross. Fig. [7| 
Lower populations of the H2(I) and H2(II) regions relative to that of the H2(III) region in 
the absence of the laser field is obvious. At the initial stages of the laser pulse, the field 
is not strong enough to alter this relative order of the populations. As the laser intensity 
increases, the outgoing population from the H2 region increases due to the single ionization 
process. The amplitudes of the oscillations of the populations of the H2(I) and H2(II) regions 
increase with increasing the intensity of the laser pulse and at some instances the populations 
of the H2(I) and H2(II) regions become more than that of the H2(III) regions. Therefore, 
the high laser field can increase population of the unstable H2(I) and H2(II) regions over 
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FIG. 7. (Color online) Time-dependent populations of different H2 regions. The left and right 
sets of panels are related to the simulations with freed and fixed nuclei, respectively. 
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that of the more stable H2(III) regions. Figure [7| also shows that when the population of 
the H~H+(I) and H+H~(II) ionic regions take over that of the covalent HH(III) regions, 
the overall population of the full H2 region decreases considerably due to the ionization. It 
worth mentioning that Dehghanian et al. [12] estimated the critical internuclear distance 
(Rc = 5a. M.) for the enhanced single ionization probability (EI) using an electrostatic model 
based on the creation of the precursor ionic state H+H^ from the covalent (natural) state 
HH occuring at the peak of the laser pulse, and showed that double-ionization probability 
has the same value as that of the single-ionization probability. The results obtained with 
fixed and freed nuclei are very similar (the left and right sets of panels of Fig. [?]) . 

QUASI-i?|+ REGIONS 

The time-dependent populations of the Q-H2''' regions are shown in Fig. [sj At the in- 
tensity of 1 X 10^^ W cm~^, there is only a small population in the Q-Hg"*" regions most 
of which is in the Q-H^+(III) region. For this intensity, the time-dependent behaviour of 
the population is the same for both cases of fixed and freed nuclei. As the laser intensity 
increases, the probability of the ionization increases and so the Q-H2''' regions become more 
populated: the populations of the Q-Il2^(I) and Q-Il2^(II) regions are increased more as 
compared to that of the Q-H2'''(III) region. In each half cycle of the laser pulse, a sharp rise 
and fall is observed for the population of the Q-H2~^(I) and the Q-H2'''(II) regions. While, the 
population variation of the Q-H2^(III) region is smoother. When the laser field is turned off, 
due to the Coulomb repulsion force between the electrons, the Q-Il2^(I) and the Q-H2^(II) 
regions are evacuated such that their population become negligible. 

H+ REGIONS 

Time-dependent populations of the singly-ionized regions H^(I) and H^(II) (Fig. [s]) are 
shown in Fig. [oj At the intensity of 1 x 10^^ W cm~^, after four optical cycles, at the 
maximum amplitude of the laser pulse (Fig. [T|, the negligible population of the single ion- 
ization region starts to increase and then becomes constant until the laser field is turned 
off. With increasing the laser intensity, populations of these regions increases drastically. 
Fig. [9] shows a correlation between populations of the H^(I) and (II) regions which can 
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FIG. 8. (Color online) The same as Fig. Q but for the Q-H^+ regions (Fig. |3| for the case of freed 
(left) and fixed (right) nuclei. 
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be regarded as a population exchange via H2 and Q-Hg intermediate regions. Comparison 
of the corresponding panels of Figs. [7]j9] shows that this population exchange occurs mainly 
through the H2 regions at low intensity and through the Q-H2''' region at high intensity. 
This population exchange becomes more evident as the intensity of the laser field increases. 
Also, the difference between the populations of the H^(I) and H^(ll) regions increases with 
the laser intensity. From Figs. |9| it is clear that at the intensity of 5 x 10^'^ W cm~^, the 
population of the singly-ionized regions is the highest compared to the other regions 
in Fig. [3] Except at 1 x 10^^ W cm~^ intensity, population of the singly- ionized regions is 
higher than that of all other regions (Fig. [s]). The populations of the singly-ionized regions 
at the 1 X 10^^ and 5 x 10^^ W cm~^ intensities show a peak in the time zones near the 
peak of the laser pulse and decreasing afterwards. This decrease is due to the outgoing of 
the population from the singly-ionized to the doubly-ionized regions ones. 



REGIONS 



Figure 10 shows the time-dependent populations of the doubly-ionized regions, (I), 
H2^(II) and H2'*'(III) (Fig. [2]). As can be seen from this figure, at the intensity of 1 x lO^'' W 
cm~^, populations of these regions are all negligible. As the laser intensity increases, these 
populations grow: the higher the intensity, the earlier the starting point of the sharp increase 
in the population. Fig. [lO] shows that the rapid rise in the population of the doubly ionized 
species occurs near the peak of the laser pulse. The populations of both doubly ionized 
H2^(I) and H2'''(II) regions are the same at low laser intensity. However, with increasing the 
laser intensity they become significantly different. Time-dependent population behaviours 
of different doubly-ionized regions are very similar for the both cases of fixed and freed nuclei 
but their values are slightly different: populations of the latter is slightly higher. 

Variation of the overall populations of the H2, H^, Q-H2^, and regions are shown 



in Fig. H As can be seen from this figure, at the lowest intensity of 1 x 10^'^ W cm~^, H2 
has the dominant population, while at the intensities of 5 x 10^^ and 1 x 10^^ W cm~^ 
has the dominant population. At the highest intensity, 5 x 10^^ W cm^^, population of the 
doubly ionized state Hj'*' is dominant. 
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FIG. 9. (Color online) The same as Fig. [7| but for the singly-ionized regions (Fig. [3| for the 
case of freed (left) and fixed (right) nuclei 
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FIG. 10. (Color online) The same as Fig. [Tj but for the doubly-ionized regions Rl+ (Fig. g for 
the case of freed (left) and fixed (right) nuclei 
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FIG. 11. (Color online) The same as Fig. [7| but for comparison of overall populations of the H2, 
H^, Q-Hg"*", and regions. The same line styles is used for different species in all panels. 



18 



POSITION EXPECTATION VALUE OF ELECTRONS 



Figure 12 shows the expectation (average) value of the position of the electrons calculated 
for two cases of fixed and freed nuclei. For both cases, this expectation value show large- 
amplitude oscillations. When the electrons under the effect of the laser field are withdrawn 
in one side from the nuclei, the attraction influence of the nuclei on the electrons decreases. 
By changing the laser field direction, the electrons come toward the nuclei and then pass 
over it in the opposite direction. This results in the oscillation of the position expectation 
value in each cycle of the laser field. At the intensity of 1 x 10^^ W cm~^, these oscillations 
loose their amplitude and become complex after turning off the laser field which ise due to 
the fact that the electronic wave function is now a linear superposition of different stationary 
states. 

At the intensity of 5 x 10^^ W cm , after the laser field termination, the resulting 
electronic wavepacket recedes from the nuclei. In this condition, some parts of the wave 
function go far away from the nuclei and result in an increase in the expectation value of 



position. Very small population transfer from to in Fig. 11 at the intensity of 
5 X 10^^ W cm~^ after 8 cycles, explains this event. At the intensity of 1 x 10^^ W cm~^, the 
population transfer are considerable such that population transfer from to starts 



earlier, about 7 cycles (Fig. 11). 

At the highest intensity of 5 x 10^^ W cm~^, a considerable population of H2 is doubly 
ionized (i.e. the electron goes out of the simulation box). In this case, with the population 
transfer from to Hg^, the reduction of the position expectation value happens. 



MOTION OF NUCLEI 

Variation of the internuclear distance due to the classical motion of the nuclei are pre- 
sented in Fig. 13 The low intensity (1 x 10^"^ W cm~^) laser field cannot overcome the 



molecule bonding energy, and thus the internuclear distance does not change considerably 
and the dissociation process does not take place. While, at higher intensities, the field is 
strong enough to break the bond of the molecule. It can be seen from Figs. [7 12 that the 



results for the freed and fixed nuclei are very similar. For the intensity of 5 x 10^^, the 
internuclear distance does not change significantly, up to cycle 6. During the cycles 6-8, 
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FIG. 12. (Color online) Variation of the position expectation value of electrons. The left and 
right sets of panels are related to the simulations with freed and fixed nuclei, respectively. 
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the internuclear distance increases slowly from ~ 3 to ~ 4.5 a.u., which leads to a small 
increase in the populations of and H^^ (thus a decrease in the population of H2) for 



the case with freed nuclei in comparison with the case with fixed nuclei (Fig. 11). For the 
intensity of 1 x 10^^ W cm~^, the internuclear distance does not change significantly up to 
the end of cycle 4. During cycles 4-8, the internuclear distance increases from 2.5 to 6 a.u. 
which corresponds to the increased ionization of the H2 and (i.e. increased population 
of H^) when the nuclei are freed (as compared to that of the fixed nuclei). This effect 
can be explained by the enhanced ionization (EI) for the internuclear critical distance (Rc) 
( \1'2\ EI])- For the intensity of 5 x 10^^ W cm~^, H2 is ionized mainly up to cycle ~ 3. At 
this intensity, the predominant population belongs to which is produced from H2 by 



double ionization and from by single ionization. Figure 13 also shows that population 
of Q-H2''' increases considerably at this intensity. The internuclear distance does not change 
significantly during cycles 0-3, but it increases rapidly from ~ 3 to ~ 8 a.u. during cycles 



3-8, as shown in Fig. 13 The increase in the internuclear distance results in a small growth 
of the Hg"*" population for the case of freed nuclei compared to that of the fixed nuclei, as it 



is clear in Fig. 11 



TOTAL NUCLEI FORCE ON ONE ELECTRON 



The total force of the nuclei on one of the electrons is extracted using 




{|fl„(i)-MP + aF^ 



This force is calculated for different laser intensities and is shown in Fig. [T4j As it is clear, 
when the laser field is on, the force is affected by the laser field oscillations so that when 
the laser pulse is turned off (at the end of cycle 8), the variations of this force becomes 
negligible. At higher intensities, the force oscillation fades out faster even before the end of 
laser field. As the laser intensity increases, the oscillation amplitude of the force increases. 
After the termination of the laser field, the distribution of electron wave packet on different 
stationary states causes a complex oscillation in the expectation value of the nuclei force on 
one electron, similar to what observed for the position expectation value of electrons. 
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FIG. 13. (Color online) Time-dependent H-H internuclear distance in the interaction of two 
electrons ID H2 system with laser field of different intensities. 

TOTAL ELECTRON FORCES ON NUCLEI 

The overall electrons exerted forces applied on the nuclei a = 1 and 2 are obtained using: 

[Ra{t) - Zi] I ij{zi,Z2,t) 12 



-'i=l 




-dzidz2 



(6) 



{[i?„(t)-^,]2 + a}3/2 

The time variation of these forces, Fee-i and Fee_2, are shown in Fig. [15} As can be seen in 
this figure, these forces are affected by the laser field oscillations and approach a constant 
value after the laser field termination. At the laser intensity of 1 x 10^^ W cm"^, since no 
dissociation occurs, these forces are the same for the two cases of moving and fixed nuclei. 
At the intensity of 5 x 10^^ W cm , dissociation of the nuclei causes the un-ionized electron 
wavepackets concentrate around each nuclei and so their effect on the other dissociated 
nucleus becomes negligible. Therefore, when the nuclei are distanced enough or when the 
electrons are ionized, the force on the nuclei is negligible. At the intensity of 5 x 10^^ W 
cm~2, the complete ionization occurs early so that the electron forces become zero at earlier 
stages of the laser pulse for both cases of the freed and fixed nuclei. The differences and 



similarities of the Fee-i and Fee-25 shown in Fig. 15, are very interesting. In contrast to the 
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FIG. 14. (Color online) Variation of the nuclei attractive force on one electron of H2 during 
and after its interaction laser field of different intensities when the nuclei are freed (left) and fixed 
(right) nuclei. 
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electrons, motions of the nuclei are distinguishable. The forces on the nuclei are symmetric 
in the absence and in the early sage of the laser field. In the higher 5 x 10^^, 1 x 10^^, and 
5 X 10^^ W cm~^ intensities, these two forces vary together to follow the variations of the 
linearly polarized laser pulse. 



TOTAL FORCE ON NUCLEI 



The total force which is exerted on the nucleus 1 is shown in Fig. 16 In comparison with 
the dimension of the simulation box, variation of the internuclear distance is very small. 
Therefore, the variations of the forces for both nuclei are very similar and we just report 
the results obtained for nucleus 1 here. In the beginning, the system rests at equilibrium 
distance. Therefore, the total exerted force on the nucleus 1 is zero. Under irradiation of 
the laser field and due to the change of the internuclear distance, the net force on nucleus 1 
does not vanish. After turning off the laser pulse, the total force on the nucleus fades out 
to zero for the freed nuclei. While, in the case of fixed nuclei, the total force never vanishes 
because of the constant repulsion force of the nuclei. For the intensity of 1 x 10^^ W cm~^, 
because of small ionization, the Coulomb attractive and repulsive forces almost cancel each 
other leading to a negligibly small total force near to zero for both cases of freed and fixed 
nuclei. 



REPULSIVE FORCE BETWEEN THE ELECTRONS 



The expectation value of the repulsive force between the two electrons is calculated by the 
first order derivative of the repulsive Coulomb potential with respect to the distance between 

dUee 

the electrons, i.e. (-- — — r). Integration is carried out over the entire computational box. 



and the results are shown in Fig. 17 Magnitude of this force reflects, inversely, the distance 
between the electrons, and therefore, can represent the amount of the correlation between 
the electrons. At the intensity of 1 x 10^'^ W cm~^, there is always a noticeable interaction 
between the electrons implying that a large part of the correlation between the two electrons 
is preserved. At higher intensities, repulsive interaction between electrons decreases rapidly, 
implying that their correlation is decreased effectively with increasing intensity of the laser 
pulse. 
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FIG. 15. (Color online) The same as Fig. 14, but for the electron forces on the nuclei 1 and 2. 
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FIG. 16. (Color online) The same as Fig. 14, but for the total force on the nuclei 1 and 2 
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FIG. 17. (Color online) The same as Fig. 14, but for the repulsive force between the two electrons 
in the simulation with freed (left) and fixed (right) nuclei. 



27 



IV. CONCLUSION 



Dynamics of the electrons and the nuclei of hydrogen molecule have been studied, based 
on the adiabatic approximation, via quantum and classical approaches, respectively by solv- 
ing time-dependent Schrodinger and Newton equations simultaneously. A one-dimensional 
model is adopted for both electrons and nuclei coordinates, and the laser-molecule interac- 
tion energy is formulated in the dipole approximation. Four different intensities have been 
used; 1 x lO^'^, 5 x 10^^, 1 x 10^^ and 5 x 10^^ W cm-^ all with 390 nm wavelength. A 
novel geometry is introduced for the simulation box used in this study. In the simulation 
boxes used in previous works, generally, the dimensions of the H2 and H2 (0^1/2 ^'^^ '^//+) 
assumed to be the same. This dimension for H2 (oJifa) should be smaller than that for 
{djj+) because the electrons of H2 must escape from a core with one positive charge, while 
the electron of must escape from a core with two positive charges. In this work, different 
values are considered for these two dimensions, which allows to reveal more details of the 
electronic wavepackct evolution and interesting properties of the system. With introduction 
of > dfj^, the regions H^(I) and H^(II) overlap in the sections named quasi-H^^ where 
both ei and 62 electrons are ionized. It is shown that there are three distinct quasi-Hg ^ 
regions. In addition, the simulation box is designed such that the components are con- 
served as long as the intensity of the laser pulse is turned on to assure that the overall 
second ionization is taken into account properly. In this research, evolution of the H2 are 
also detailed by dividing it into four sub-regions related to the ionic H+H^ and covalent 
(natural) HH states. It is shown that when the population of the ionic regions H~H"'" (I) 
and H+H" (II) becomes higher than that of the covalent regions HH(III), ionization causes a 
considerable reduction in the population of the H2 region. The time-dependent populations 
of different regions are calculated at different intensities and analyzed comparatively. It is 
shown that at the lowest intensity of 1 x 10^^ W cm~^, H2 has the dominant population, 
while at the intensities of 5 x lO'^^ and 1 x 10^^ W cm"^, has the dominant population. 
At the highest intensity, 5 x 10^^ W cm~^, the population of the double ionization region, 
H2''', is dominant. It is shown that the rapid increase of the double ionization population 
occurs near the peak of the laser pulse. The double ionization population of both H2'''(I) and 
H2''"(II) regions are the same at the lowest intensity examined here. With increasing the laser 
intensity, the population of the H2''"(I) and H2'''(II) regions becomes significantly different. 

28 



The H2 system does not proceed to dissociation at the lowest intensity of 1 x 10^ W cm~^, 
but it is dissociated at higher intensities. Also, the effect of the internuclear distance and 
motion of nuclei on the enhanced ionization is discussed. Finally, different time-dependent 
properties of the system were calculated and analyzed based on the characteristics of the 
laser pulse and variation of the populations of different regions. These properties include 
the internuclear distance, the total force of the nuclei on one of the electrons, the total 
electron forces on nuclei, the total force exerted on the nuclei, and repulsive force between 
the electrons. 
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